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We report on results from three spectrum calculations with staggered quarks: 1) a quenched calculation with 
the standard action for the gluons and quarks; 2) a quenched calculation with improved actions for both the 
gluons and quarks; and 3) a calculation with two flavors of dynamical quarks using the standard actions for the 
gluons and quarks. 



1. INTRODUCTION 

A controlled calculation of the masses of the 
light hadrons has been one of the major goals of 
lattice QCDjiJ. It has become increasingly clear 
in recent years that in order to carry out such a 
calculation one must perform high statistics simu- 
lations using a range of lattice volumes in order to 
control finite size effects, a range of quark masses 
in order to extrapolate to the chiral limit, and 
a range of lattice spacings in order to extrapo- 
late to the continuum limit. We have carried out 
such a series of simulations in the quenched ap- 
proximation with staggered quarksj|]. In Sec. 2 
we report on new results from this study with 
particular emphasis on the extrapolation to the 
chiral limit. 

Improved actions have the potential to signif- 
icantly increase our ability to perform accurate 
spectrum calculations |J. In Sec. 3 we report 
on exploratory quenched calculations using im- 



proved actions for the gluons and the staggered 
quarks. Finally, in Sec. 4 we show preliminary re- 
sults from a series of full QCD calculations with 
two flavors of staggered quarks at gauge coupling 
6/ ff 2 = 5.50. 



2. QUENCHED SPECTRUM 
STAGGERED QUARKS 



WITH 
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Over the past several years we have carried 
out a series of quenched spectrum calculations 
with staggered quarks. Calculations were per- 
formed with four values of the gauge coupling, 
6/g 2 = 5.54, 5.70, 5.85, and 6.15. Five differ- 
ent lattice volumes were used at 5.70 and three 
at 5.85 in order to study finite volume effects. 
Five quark masses were used at each coupling. 
For 6/g 2 = 5.70 and 5.85, we used am q — 0.16, 
0.08, 0.04, 0.02 and 0.01. At 6.15 the masses were 
halved, and at 5.54 they were doubled. For each 
coupling, lattice volume and quark mass we per- 
formed correlated fits to determine the masses. 
An Edinburgh plot is shown in Fig. 1. 



2 



1.6 



* 5 



© 

to 



1.4 



6/ 



1.2 



* 5.54 16 J 

« 5.70 20 3 

O 5.70 24 3 

X 5. 85 20 3 

□ 5.85 24 3 

o 6.15 32 3 



0.0 



0.5 

m„/ m. 



1.0 



Figure 1. Edinburgh plot for the quenched stag 
gered spectrum. 
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Figure 2. /m p as a function of lattice spacing 
for fixed m^/rrip. 



Hadron masses calculated on the same set of 
gauge configurations, but with different quark 
masses are, of course, correlated. So, we per- 
formed correlated fits for each coupling and lat- 
tice volume to extrapolate to the chiral limit. We 
used a wide variety of fitting functions with terms 
inspired by quenched chiral perturbation theory. 
The first conclusion was that simple linear or 
quadratic fits to the nucleon and rho masses are 
inconsistent with our data. Good fits are possible. 
For example, the function M + am q + bm q ^ 2 + cm q 
gives reasonable fits for both the nucleon and 
rho masses. We also obtain good fits for the 
nucleon with functions containing square roots 
or logarithms of the quark masses suggested by 
quenched chiral perturbation theory. Once such 
a fit is in hand, we can investigate the dependence 
of tojv /Trip as a function of lattice spacing for fixed 
m^/rrip. Some typical plots are shown in Fig. 2 
using the above fitting form. These fits allow us 
to extrapolate to the continuum limit for fixed 
values of m^/rrip. They are also proving very use- 
ful in comparing results obtained with improved 
actions to our high statistics results with the stan- 
dard staggered action. 

In Fig. 3 we plot m^/mp as a function of 
arrip(m q = 0) with the ratio m^/rrip held fixed 



at its physical value. Fits to the rho and nucleon 
used in this plot were again made with the form 
M + am q + brr^J 2 + cm 2 ,. The solid curve is a 
fit of the form c + c 2 a 2 to data for all four cou- 
plings, and the dotted line is a fit of the same form 
to data for the three weakest couplings. Here a 
is the lattice spacing, and Cj are fitted parame- 
ters. We are unable to uniquely determine the 
optimum fitting function for the chiral extrapo- 
lation from our data. In particular, nucleon fits 
with m)/ 2 terms have confidence levels similar to 

1 /2 

those shown in Figs. 2 and 3. An m q terms tend 
to pull the value of rriN down in the chiral limit, 
thereby reducing mjv/mp. At present the uncer- 
tainty in the choice of fits is the largest source of 
error in our determination of rrifq /m p . 

3. QUENCHED SPECTRUM WITH IM- 
PROVED STAGGERED QUARKS 

Up to now, most of the effort aimed at im- 
proving the quark action has focused on Wilson 
quarks; however, some time ago Naik proposed 
adding third nearest neighbor hopping terms to 
the staggered quark action, adjusted to remove 
the order a 2 errors in the lattice quark propa- 
gatorfl. We are experimenting with a tadpole 
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Figure 3. jm p as a function of am p (m q — 0) 
with the ratio m^/rrip held fixed at its physical 
value. 



Figure 4. The Edinburgh plot for the quenched 
spectrum with improved gluon action at [3 p i = 
7.4. 



improved version of this action to carry out a 
quenched spectrum calculation with gauge config- 
urations generated using a tadpole-Symanzik im- 
proved gauge action^. We have generated con- 
figurations at three different couplings, f3 p i = 6.8, 
7.1 and 7.4, on lattices ranging in size from 
8 3 x 16 to 16 3 x 32. We have measured the 
spectrum on these lattices with both the tadpole 
improved Naik action and the standard Kogut- 
Susskind action. In Fig. 4 we show the Edin- 
burgh plot for the Naik (crosses) and standard 
staggered (diamonds) quark actions on 16 3 x 32 
lattices at (3 p i — 7.4. For comparison, we also 
show quenched spectrum results with the stan- 
dard gauge and staggered fermion actions (fancy 
plusses) at 6/g 2 — 5.54. This coupling corre- 
sponds to approximately the same lattice spacing, 
as measured by the rho mass, as the improved ac- 
tion at P p i = 7.4. One sees from Fig. 4 that the 
improved gluon action does lead to lower values 
of the nucleon to rho mass ratio at small values of 
m^/rrip, than the standard gluonic action. How- 
ever, the Naik term has little effect at the values 
of the coupling constants and quark masses that 
we have investigated to date. 

An important test of improved staggered quark 



actions is the extent to which they restore fla- 
vor symmetry. Our preliminary results indicate 
that the mass difference between the ir and 7T2 
does decrease with the improved gluon action, but 
that the Naik term has little impact. Although 
our tests of the tadpole improved Naik action 
have not been very encouraging at the parameters 
studied to date, it should be noted that this ac- 
tion has the advantage of significantly hastening 
the rate at which thermodynamic quantities ap- 
proach the continuum limit for free quarks. One 
might therefore hope that it would lead to a sig- 
nificant improvement in high temperature QCD 
calculations, and the Bielefeld group has shown 
that this is the case||. It would therefore seem 
that this action should be studied in greater de- 
tail. 

4. SPECTRUM WITH TWO FLAVORS 
OF STAGGERED QUARKS AT 6/g 2 = 
5.50 

As part of our study of the decay constants of 
heavy-light mesons, we are generating a set of lat- 
tices with two flavors of Kogut-Susskind quarks at 
gauge coupling 6/g 2 — 5.50. We are using quark 
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masses m q = 0.10, 0.05, 0.025 and 0.0125. All lat- 
tices have a time dimension iV t = 64. We use spa- 
tial volumes 20 3 for the two lighter masses, and 
24 3 for the two heavier ones. (The use of larger 
spatial volumes for the heavier quark masses was 
dictated by the peculiarities of the CM5 on which 
the heavier quark mass calculations were per- 
formed). We have calculated the spectrum us- 
ing the first 1,000 equilibrated molecular dynam- 
ics time units for each quark mass. The Edin- 
burgh plot is shown in Fig. 5. We have included 
quenched data for gauge couplings 6/.g 2 = 5.70, 
5.85 and 6.15 for comparison. The data indicates 
a somewhat faster fall off of tojv/to p than the 
quenched data at 6/<? 2 = 5.70, which has a com- 
parable lattice spacing. This is to be expected 
because lighter dynamical quarks will lead to a 
greater renormalization of the coupling. 



trajectories at each mass. 
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Figure 5. The Edinburgh plot for full QCD with 
two flavors of staggered quarks at 6/g 2 = 5.50. 

In Fig. 6 we plot the masses of the tt and the 
7T2 as a function of am p to illustrate the extent 
of flavor symmetry breaking in this data set. We 
also include very preliminary results from a set 
of runs at 6/g 2 = 5.60 on 24 3 x 64 lattices with 
quark masses 0.01, 0.02, 0.04 and 0.08. These 
results are from approximately 600 equilibrated 
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Figure 6. m v and m I2 as a function of am p for full 



QCD at 6/g 2 = 5.50 and 5.60. 
are quite preliminary. 



The 5.60 results 



The calculations described above were car- 
ried out on Intel Paragons at Indiana University, 
ORNL and SDSC, the TMC CM5 at NCSA, the 
Cray T3D at PSC, and DEC Alphas at Indiana 
University and PSC. This research was supported 
by the United States NSF and DOE. 
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